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Abstract 

The contribution Ae of extremely small motions of domain walls to small-signal permittivity of 
a multidomain ferroelectric sample has been a research issue for many years. In ferroelastic ferro- 
electrics such motions contribute also to their piezoelectric (by Ad) and elastic (by As) properties. 
Data about their simultaneous existence are scarce but those available point to mutual propor- 
tionality of Ae, Ad and As, as expected. To understand the magnitude of extrinsic contributions, 
the origin of the restoring force acting on domain walls must be understood. In the present con- 
tribution the theory has been developed based on the model of a plate-like sample in which the 
ferroelectric- ferroelastic bulk is provided with a nonferroic surface layer. Motion of domain walls 
in the bulk results in a change of electric and elastic energy both in the bulk and in the layer, which 
provides the source of restoring force. This makes it possible to determine all mentioned extrinsic 
contributions. We discuss the applicability of the model to available data for single crystals and 
also for ceramic grains. 
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I. INTRODUCTION 



The problem of extrinsic (domain wall) contributions has been addressed by many au- 
thors, both experimentally and theoretically. In the prevailing number of cases, only extrin- 
sic permittivity has been studied. For piezoelectric ceramics, Arlt et al.|l| were the first to 
address the problem of wall contributions to all involved properties: permittivity e, elastic 
compliances s and piezoelectric coefficients d. In this and related papers, the existence of the 
restoring force is assumed and its origin was not specified. Later, Arlt and Pertsev|2| offered 
a more involved approach: when domain wall in a ceramic grain is shifted, uncompensated 
bound charge appears on the grain boundary, producing electric field. Simultaneously, if 
the involved domains are ferroelastic, a wall shift results in mechanical stress in surrounding 
grains. Thus the shift is accompanied by the increase of both electric and elastic ener- 
gies, leading to restoring forces. Results of these theories have been successfully related to 
experimental data on all three mentioned contributions: Ae, As anc 

nnn 

While there exist a number of such data for ceramic materials[ll, 13|, information for sin- 
gle crystals is rather scarce. Understandably, for nonferroelastic ferroelectrics such as TGS 
only data on Ae are available; a] on the other hand for crystals which are both ferroelectric 
and ferroelastic with more than two domain states, dense domain systems are rather chaot- 
ical, difficult to approach theoretically. In the present paper, we have in mind ferroelectric 
and ferroelastic crystals with only two domain states. In particular, crystals belonging to 
the KH2PO4 family belong to this category and have been intensively studied. Nakamura et 
al.le,!?] determined Ae^s for KH2PO4, CSH2PO4 and CSH2ASO4. For KH2PO4, Nakamura 



and Kuramoto[8[ proved the existence of both Ae^s and Asqq while Ad^Q was measured for 
RbH,PO J. Fo. the =a„.e .naterial, all ,uantM« A<*. and have been n.easu.ed 

by Stula et al.[10[ It was found that all these contributions are mutually proportional when 
measured as a function of temperature, in the temperature interval between Tq and Tc-35 
K. Several other ferroics for which our approach may be applicable will be mentioned at the 
end of this paper. 



In single crystals, the origin of the restoring force is usua^ 



pinning on crystal lattice defects. In our recent papers ll|, ll2(] we introduced the model of 



ly connected to domain wall 



a passive surface layer to calculate the restoring force for domain walls and the resulting 
extrinsic permittivity and piezoelectric coefficient. In fact, the infiuence of a surface layer 
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FIG. 1: Geometry of the model in x-z plane 



on the properties of a ferroelectric sample was discussed repeatedly several decades ago. In 
particular, in connection with the sidewise motion of domain walls in BaTiOs, thickness 
dependence of the coercive field, asymmetry of a hysteresis loop or the problem of energies 
of critical nuclei, theoretically impossibly high. 

In the present paper, we return to this approach. However, in contrast to previous 
calculations, we offer a more involved model. The shift of a domain wall induced by the 
application of electric field or elastic stress results in the increase of both electric and elastic 
energies. In the following, these are explicitly calculated which makes it possible to determine 
all extrinsic coefficients Ae33, Asqq and Ad^Q. Their dependence on the sample properties 
will be discussed. 



II. DESCRIPTION OF THE MODEL 



We consider a plate-like sample elecroded sample of infinite area with major surfaces 
perpendicular to the ferroelectric axis z. The material is both ferroelectric and feroelastic; 
domains with antiparallel polarisation differ in the sign of spontaneous shear. However, we 
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FIG. 2: Geometry of the model in x-y plane 
shall approximate the material in the ferroelectric phase by the equation of state 

Dx = S^SaExi 

= EqEcE^ + Po, 

neglecting the intrinsic piezoelectricity. As we shall see, this does not supress the existence 
the extrinsic piezoelectricity, which is one of the aims of our calculations. In the preceding 
equation, we also neglect nonlinear terms; for 2nd order phase transitions this limits the 
validity of our calculations to the temperature region not very close bellow the temperature 
Tq- Domain walls are assumed to have surface energy density a^, and zero thickness. 

For simplicity, we approximate elastically anisotropic material of the sample by the elas- 
tically isotropic one. Neglecting again the intrinsic piezoelectricity, its mechanical properties 
are described by equations for stress tensor components 

where /i^ and fif are Lame coefficients of the passive layers and the bulk respectively; Cij 
is the strain tensor, eo,jj is the spontaneous strain tensor of the central ferroelastic part. 
We suppose that the only nonzero components of the spontaneous strain tensor are eo,i2 = 
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eo,2i = ±eo in the a_|_ resp. a_ domain, see Fig. [21 We introduce the asymmetry factor 



We neglect thermal interactions and suppose that the sample is thermally isolated. To 
keep the constant voltage V on the sample, the electrodes should be connected to external 
electrical source. In the same way, to keep constant external stress Text, 12 = Text, the sample 
should be deformed by external mechanical force. The infinitesimal work of these external 
sources should be taken into consideration when discussing the variations of the energy of 
the isolated system sample + sources. 

III. HELMHOLTZ FREE ENERGY 

In what follows we consider three contributions to Helmholtz free energy, calculated per 
unit area of the plate-like sample (in the x-y plane): 

The energy of domain walls per unit area of the sample 

the electric field energy per unit area of the sample 



where the integration is taken over parallelepiped x G (0, 2a), z G (— t/2), y G 
(0, 1 m^/2a), and energy of elastic deformations per unit area of the sample: 



where the integration is taken over the same region, in the first case for constant external 
stress Text in the plane x = 2a, in the second one for constant displacement u in the plane 
X = 2a. In both cases, the displacement for x = is chosen to be zero and the boundaries 
of the sample in x-y plane are free of stress. 

To find the Uei and Udei, we have calculated electric potential and mechanical displacement 
inside the sample by Fourier method. We present here only the relative simple results for a 
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resp. 
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"dense pattern approximation" , that is for d, h ^ a: 

dhA^P^+ V^elsdSc 



UeiiV, A) 



f^def(Text, A) 



Udei{u, A) 



2eo {Sd h + Ecd) 



2{iidd + fifh) 2 



fidd + fifh 



Infinitesimal work of the electric source at constant voltage V per unit area of the sample 
plate is 

6W^i = V 6ao [J ■ m-2] , 

where ctq is constant Fourier component of the free charge density on the positive electrode, 
calculated as 

PqA + EqE^ V/h 



ao{V, A)- M • 

1 + Bed/ [Edh) 

Infinitesimal work of the mechanical source, deforming the parallelepiped is 



(1) 



It is easy to prove that the equilibrium domain structure for V = Text = resp. 
u{x = 2a) = is symmetric (i.e. A=0), with domain width 



Oeq = V 3.68 ha^ 



4:6, 



2 f^df^f 



p2 



fid + /i/ eo {Sd + y/Sa Sc) 



(2) 



It can be shown (see e.g. ll|]) that within a large interval of the applied voltage the average 
width 

a = (a+ + a ) /2 

remains constant. This is why can be also considered as constant. 



IV. EXTRINSIC CONTRIBUTIONS Aegg, Asee AND Ad, 
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We calculate the equilibrium effective egf of the sample from the relations 



Df = aoiV, A) = ef, E^' = el^ 



V 



(3) 



We keep Text = 0, V^=constant and we take into account that variation of "Helmholtz free 
energy of the sample + the work of electric source" is zero in equilibrium: 



dU,i{V, A) 



5A 



<9t/def(rext 



0, A) 



5 A - 5Wei = 0. 



dA dA 
Solving this standard problem, we get A = A{V) and from the Eqs. (Q, the effective 



£33. For fid = fJ'f = fJ', Sd = £c = £z obtain the relatively simple result: 



-eff 



£z + £z 



h 
d 



P^h + Ael EoEz /it 



(4) 



The effective elastic compliance of the sample is 



eff _ 4 eff 



o eff 



U 



1212 



Text 2 a Text 

We put V = (shorted sample), apply constant external shear stress Text and postulate, 
that variation of "Helmholtz free energy of the sample + the work of mechanical source" is 
zero in equilibrium 

oA oA 

Solving this problem, we get u = M(rext) and we get for fid = fij = jU, Ed = £c = £z- 



.eff 
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H jj, d 



4 el Eq e. 



To find the effective piezoelectric coefficient of the sample 



(5) 
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r)eff 
^3 



Text, 6 Text, 6 

we put V = 0, apply Text and solve the problem 

a^ei(0, A) , dU,,,iu, A) , dU,,,iu, A) ^ 

From here we obtain A = y4(rcxt)- Inserting this result into Eq. ([T]) we obtain cxq = o"o(T( 
For iJid = Hf = lJ^,Ed = Ec = Ez^e get finally for the effective piezoelectric coefficient 

reff _ ^ 2eo Pq Eq E^ 

d' [p,^ + AelEoEzfi_ 
The same result we get for the inverse piezoelectric effect. 



cxt J 



(6) 
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TABLE I: Numerical estimate of Ae33, Asqq and Adsg for different values of surface layer thickness. 
The following numerical constants have been used: Pq = 4-10~^Cm~^, = 100, eg = 0.015, 
fi = e-lO^Pa, 0-^ = 5-10-3 J m-2, h = 5-10-^ m. 





d = 0.5-10-6 m 


ci = 2.5- 10-6 m 


d = 5-10-6m 


Aeas [1] 


12 000 


2 500 


1200 


Ad36 [Cm-2] 


4.1-10-8 


8.3-10-9 


4.1-10-9 


As66 [Pa-^] 


1.6-10-8 


3.2-10-9 


1.6-10-9 



V. DISCUSSION 



The above calculation leads to explicite formulae (jl]) to ([6]) for extrinsic permittivity, 
extrinsic elastic compliance and extrinsic piezoelectric coefficient, resp. Numerical values 
for all involved material coefficients are available for single crystals of RbH2P04: Pq = 
4-10-^Cm-^, 6z = 100, Co = 0.015 and fi = 6 -109 Pa. To obtain numerical estimates for 
particular samples we put = S-IO-^ Jm-^ and h = 5-10-^m and choose three values of 
surface layer thickness, namely d = O.S-lO-^m, d = 2.5-10-6m and d = S-lO-^m. Table [T] 
shows resulting values of all three extrinsic variables. These numbers appear very realistic 
and confirm the applicability of the model presented in this paper. 

It is appropriate to pay some attention to the fact that also the formula ([2]) gives a rea- 
sonable numerical estimation for the width of equilibrium domain pattern. With numerical 
values specified at Tab. [T]we obtain Oeq = 1 /im. 

In the approach analyzed above, the source of the restoring force acting on domain walls 
is the interaction of ferroic sample with a passive surface layer. Very often, the origin of 
restoring forces is connected with domain wall pinning to defects. Understandably, the 
latter mechanism cannot be excluded for ferroics of any chemical composition. On the other 
hand, passive surface layers can be formed during sample preparations and in particular 
for water-soluble crystals their appearance is a very likely: samples are polished in water- 
containing media, the procedure obviously leading to the presence of a passive surface layer. 
The example analyzed numerically above, crystals of RbH2P04, falls into this category. 
However, extrinsic properties of a number of other crystals have been studied. Thus, e.g.. 
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for LiTlC4H406 ■ H2O (species 222-Peds-2y) [13|] very strong and nonhysteretic dependence 
of sf^ as a function of applied field E was measured. This is possible to explain by a strong 
contribution of domain walls with a pronounced restoring force. The above model would 



lead to such behavior. Similarly, large extrinsic contributions to sn have been measured 14] 
for (NH4)4LiH3(S04)4, species 4-eds-2. 
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